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Abstract 

To investigate relations between long-range antiferromagnetic (AF) order, 

superconductivity and two particle triplet collective excitations we consider a 

modified two dimensional t—J model at doping close to half filling. The model 

includes additional hopping t" and nearest sites Coulomb repulsion V . The 

additional parameters allow us to control closeness of the system to the AF 

instability. We demonstrate the possibility of co-existence of long-range AF 

order and d-g-wave superconductivity. In the phase with long-range AF order 

we find, analytically, superconducting gaps and spin wave renormalization. 

We demonstrate that at approaching the point of the AF instability the spin 

triplet collective excitation arises with energy below the superconducting gap. 
PACS: 71.27.+a, 74.20.Hi, 75.50.Ee 
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It is widely accepted now that superconductivity of cuprates is closely related to their 
unusual magnetic properties, and it is increasingly clear that magnetic pairing is the most 
realistic mechanism of cuprate superconductivity. However the mechanism of pairing as well 
as other unusual properties are far from completely understood. The problem has been 
attacked along several directions. First we have to mention the empirical or semi-empirical 
approach which allows one to relate different characteristics measured experimentally. This 
approach is to a large extent based on the Hubbard model. For a review see articleB. In the 
low energy limit, the Hubbard model can be reduced to the t— J model. Another approach to 
cuprates is based on numerical studies of the t — J model (see review!). Our studies are also 
based on this model. We used the ordered Neel state at zero doping as a starting point to 
develop spin- wave theory of pairingll. The method we used was not fully satisfactory, since it 
violated spin-rotational symmetry, nevertheless it allowed us to calculate from first principles 
all the most important properties including critical temperature, spin-wave pseudogap and 
low energy spin triplet excitations^. 

A sharp collective mode with very low energy has been revealed in YBCO in spin po- 
larized inelastic neutron scatteringjUl. A number of theoretical explanations have been 
suggested for this effectiS. All of those explanations are based on the idea that the system 
is close to AF instability. However, all known explanations use some uncontrolled approxi- 
mations and assumptions. 

In the present work we investigate a close to half filling regime for the 2D t — J model 
where it can be solved analytically without any uncontrolled approximations. It can be 
done for the region of parameters where long-range AF order is preserved under doping. 
We analyze the superconducting pairing in this regime and consider spin triplet collective 
excitation. It is demonstrated that the narrow collective mode arises at approaching the 
point of AF instability and the energy of the mode is below the superconducting gap. 

Let us consider a,t — J — J" — V model defined by the Hamiltonian 



H — t ^ c|g.CjCT t ^ c|^^Cj2o- -|- ^ 

(ij)a- (ij2)o- (ij) 



J ( SjSj — -niUj ) -|- VniTij 



(1) 



is the creation operator of an electron with spin cr (cr =1, |) at site i of the two-dimensional 
square lattice. The cj^ operators act in the Hilbert space with no double electron occupancy. 
The (ij) represents nearest neighbor sites, and {ij2) represents next next nearest sites. The 
spin operator is Sj = | Z]a,/3 c|a'^a/3Ci/3, and the number density operator is rii = J^aclaCia- 
In addition to the minimal t — J model (see Ref.i) we have introduced additional next next 
nearest hopping t", and Coulomb repulsion V at nearest sites. Note that we do not introduce 
next nearest neighbor hopping t' (diagonal) because we do not need it for the purposes of 
this study. 

In the present work we consider the perturbation theory limit t, t" ^ J. It is well 
known that the t — J model with small hopping has phase separation at doping. The reason 
for this is very simple: attraction between holes at nearest sites due to the reduction in 
number of missing AF links. The value of this attraction immediately follows from eq.(0): 
U ~ J(SjSj — 1/4) ^ —0.58 J. This effect is very simple and is not of interest to us. To 
overcome it we choose a Coulomb repulsion V = 0.58 J. So the only purpose of introducing 
the Coulomb repulsion in the model is to eliminate short range hole-hole attraction. 

We consider phase with long range AF order. It is well known (see e.g. Ref.i) that minima 
of the single hole dispersion are at the faces of the magnetic Brillouin zone (±7r/2, ih7r/2), 
and the dispersion is (we take energy at the minimum as a reference point) 

ek = f3ill + f32{l^)^ (2) 
Pi = 3.26t^/J + 8t", P2 = l.2lt^/J + 8t". 

— cosfcy). Calculation of the t" contribution in the 
dispersion (0) is straightforward because it is hopping within the same magnetic sublattice. 
Calculation of the t-contribution is more complicated, and we use the results of a series 
expansion^. Note that a selfconsistent Born approximation0 gives substantially different 
coefficients {(3i = 1.86t^/J, (32 = 0.26t^/J). The reason is that there are several important 
diagrams at small t which are missed in this approximation. We will consider the case of 
very small doping, 5 -C 1, with respect to half filling (total filling is 1 — 5). In this case 
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all holes are concentrated in small pockets around the points ko = (±7r/2, ±7r/2). Single 
hole dispersion (g) can be expanded near each of these points = Pipf/2 + (32p\/'2,, where 
p is deviation from the center of the pocket: p = k — ko, pi is orthogonal to the face of 
the magnetic Brillouin zone, and p2 is parallel to the face. The Fermi energy and Fermi 
momentum for the holes equal ep ~ |7r(/5i/?2)^''^'^, Vf ~ ^JPifP2F ~ (tt^)-*^/^. 

Spin-wave excitations on an AF background are usual spin waves with dispersion cUq = 
2 J^l - 72 ^ V2Jg, at g << 1, see Ref.li=l for review. The hole-spin- wave interaction is well 
known (see, e.g. Ref.0) 



Hh^sw = J2 ^k,q ( /^k+qi^kTttq + ^k+qT^ki/3q + H.C. ) , (3) 
k,q ^ ^ 

5fk,q = 4ty2(7kt/q + 7k+qV;j), 

where hl_^ = Ck,-o- is the hole creation operator, a], and M are the spin wave creation 



operators for = =f1, and f/q = \J^- + \ and = —sign{'^ci)\J-^ — \ are parameters of 
the Bogoliubov transformation diagonalizing spin-wave Hamiltonian, see To describe 

renormalization of the spin wave under doping, it is convenient to introduce the set of Green's 
functionJii 

l)aa(t,q) = -^(r[aq(t)at(0)]), (4) 

D„^(t,q) = -z(r[aq(t)/5_q(0)]), 
D^„(t,q) = -«(r[/?lq(t)at(0)]), 

In the present work we consider only the long-range dynamics: q ^ k pp <^ 1. In 
this limit all possible polarization operators coincide! Paa{oj^c\) = Pai3{uj,q} = P^a(u;,q) = 
Pj3f3{u!,q) = n(ct;,q). For stability of the system the condition (Stoner criterion) 

ujg + 2n(0, q) > (5) 

must be fulfilledlli. Otherwise the Green's functions (^) would possess poles with imaginary 
uj. Considering holes as a normal Fermi liquid one can easily calculate the polarization 



operator at g ^ pp: 11(0, q) ~ — 4t^-\/2g/7rv^/3i/52, Ref.0. Relatively weak pairing, which 
we consider below, does not influence this result. Then the condition of stability can be 
rewritten as 



4K^2 > ^. (6) 



ttJ 

The dispersion (|^) at t" = violates this condition. This is a well known statement about 
the instability of long-range AF order under doping in a pure t — J model. To provide 
stability we have to choose t" > = 0.0635t^/J. We want to have the system in the phase 
with long range AF order, but close to the point of instability, and therefore we will consider 
t" only slightly above the critical value. As a measure of this closeness it is convenient to 
introduce the parameter rj 

The criterion (|^) is proportional to this parameter. 

Superconducting pairing mediated by the spin-wave exchange can be found analytically!. 
Only the pairing within one pocket is important. There is no solution for the gap without 
nodes (s-wave), and there is an infinite number of solutions with the nodes. The pairing is 
maximum for the case of a single node line in the pocket, and the gap at the Fermi surface 
(ep = (3ipI/2 + f32pl/2) is of the form 

A(0) = Aosin0, (8) 

Ao = Ce^e-i/^ 

where 



sm = ^ (9) 



and C ~ 1 is some constant. Energy spectrum and Bogoliubov parameters are given by the 
usual BCS formulas 



E^ = ^{e,-eFf + Al (10) 

Substituting values /5i = "i.llt^ j J and /52 = l.lSt^ / J corresponding to the critical point, we 
find g ~ 0.25, and hence 

Ao ^ epe"^ ^ 0.07^. (11) 

The eqs.(^,(P), and (|r^) describe pairing within a single pocket. There are effectively two 
pockets in the magnetic Brillouin zone. Taking symmetric and antisymmetric combinations 
between pockets we get d- and g-wave pairings respectively (see for details). The gaps, 
as well as the critical temperatures for d- and g-wave, are practically the same. Note that we 
discuss the situation without short range hole-hole interaction, which has been eliminated 
by adjusting the nearest sites hole-hole Coulomb repulsion {V ~ 0.58 J, see above). If 
one destroys this adjustment, the degeneracy between d- and g- waves is also destroyed. 
The crucial point is that the g-wave is not sensitive at all to the interaction at nearest sites, 
whereas the d-wave is very sensitive. Therefore at < 0.58 J the d-wave pairing is enhanced, 
while on the contrary, at larger Coulomb repulsion V > 0.58J the d-wave is suppressed and 
disappears very fast. 

Now we can discuss the spectrum of the spin-triplet excitations. The Spin wave polar- 
ization operator due to the mobile holes is of the form (see e.g. Ref.S) 

n(^, q) = 2^ ^koq^ , F ^ ( Wk^k+q + Mk^^k^k+q^^k+q ) • (12) 

k,ko ^ - l^k + ^k+qj ^ ^ 

It accounts for both normal and anomalous fermionic Green's functions. Eq. ([12D includes 
summation over pockets ko = (vr/2, ±7r/2). In these pockets the vertex is g]^o,^ ~ 2^/'^)f:(g^. ± 
qy)/y/q- Let us consider the case of very small momenta and frequencies: Vpq < Aq, and 
uj < Aq. In this limit one can put g = in eq. (|l^) everywhere except at the vertex and 
therefore the polarization operator can be evaluated analytically 
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Note that the the imaginary part is nonzero even at u < 2Aq because the gap (P) has a hne 
of nodes. Any of the Green's functions (^) have a denominator cj^ — cj^ — 2a;qn(ti;, q), see 
e.g. Refs.00. The zero of this denominator gives the energy and width of the spin-triplet 
collective excitation. Using eqs.(|T3|) and we find 

Oq = ?7^q, (14) 
= 8Ao"-- 

In essence this is the renormalized spin-wave, but its energy is much smaller than the energy 
of the bare spin- wave, Oq/wq = 1] <^ 1. This collective excitation exists as a narrow peak 
only at very small q when vrc^q/SAo < 1. At higher q the width is larger than its frequency 
because of decay to particle hole excitations. 

The last question we want to discuss is how close to the point of AF instability can 
we approach using the present description, or in other words how small can rj be? The 
superconducting pairing itself is not sensitive to the AF instability and it survives at the 
transition to the disordered phase. The reason for this is clear: even in the disordered phase 
the magnetic correlation length C,m is much larger than the typical distances r ~ 1/pf ~ 
l/(7r5)^/^ which are important for pairing. However the spin-triplet spectrum (Q) is valid 
only if 1] is not too small. The sublattice staggered magnetization m is renormalized because 
of spin-wave spectrum renormalization, see Ref.i 

J uj (277)2 V ) 

The integral converges at the momenta q where the collective excitation exists as a narrow 
mode: q < Aq/J. Our consideration is valid if the renormalization of staggered magneti- 



zation is small: 6m ^ 0.5. Together with eq. (|15]) this gives the limit for validity of our 
approach. 

V > ^o/J (16) 



At Aq > T]/ J > equation ([1^) for the spectrum is not valid, but nevertheless there is some 



gapless spin-triplet excitation. The point 77 = indicates transition to the spin hquid phase 
and a spin-wave gap is opened here. 

In conclusion we have considered a close to half filling t — J — J" — V model with Coulomb 
repulsion y = 0.58J adjusted to ehminate short range hole-hole interaction. We restrict our 
consideration by the case of small t,t <^ J, and small doping 5 <^1. It is demonstrated that 
at t" > 0.0635t^/J the Neel order is preserved under the doping, and at t" < 0.0635t^/J 
the order is destroyed and the system undergoes transition to the spin liquid phase. In both 
phases there is d-g-wave superconducting pairing mediated by spin-wave excitations. In the 
Neel state we found, analytically, collective spin triplet excitation. It exists as a narrow 
mode only at very small momenta and its energy is substantially below the energy of the 
bare spin wave. 

We wish to thank M. Kuchiev for stimulating discussions. 
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